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Abstract. In this paper we consider the motion of a Brownian particle in an 
inhomogeneous environment such that the motion can be described by the equation 
yielding 1// spectrum in a broad range of frequencies. The inhomogeneous 
environment can be a result, for example, of a linear potential affecting the Brownian 
particle together with the medium where steady state heat transfer is present due 
to the difference of temperatures at the ends of the medium. The correlation of 
collisions between the Brownian particle and the surrounding molecules can lead 
to the situation where the finite correlation time becomes important, thus we have 
investigated the effect of colored noise in our model. Existence of colored noise leads to 
the additional restriction of the diffusion and exponential cut-off of the distribution of 
particle positions. Narrower power law part in the distribution of the particle positions 
results in the narrower range of frequencies where the spectrum has power law behavior. 
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1. Introduction 

One of the characteristics of the noise is the power spectral density (PSD). The white 
noise has a frequency-independent PSD, whereas the PSD of a colored noise depends 
on the frequency, the characteristic behavior of the PSD is referred to as a “color” of 
the noise. There are various applications where the noise in the physical system under 
investigation has a non-trivial spatio-temporal structure and where it is not realistic to 
model it is as a white noise process. For example, for a Brownian particle the driving 
noise is actually colored, i.e. it has a characteristic non-zero correlation time r on a short 
time-scale of the order of tens of nanoseconds jl]. Noise color arises due entrainment 
of fluid around diffusing particle. The particle accelerates the entrained fluid and this 
acceleration depends on the past motion of the particle and introduce an inertial memory 
effect [2]. 

Theoretical models suggests that variety of systems affected by a colored noise 
instead of a white noise exhibit new interesting properties. For example: The 
correlations in colored noise are found to be able to enhance or suppress the growth 
rate of amplihcation above or below a critical detuning in the collective scattering 
of light from a laser with a colored noise in ultracold and collisionless atomic gas 
[3]. Investigation of the colored-noise effect on nonequilibrium phase transitions shows 
reentrant transitions from ordered into the disordered phase as the correlation time 
and the coupling strength increase [1]. The color and coupling induced disorders are 
pure colored-noise effects because of the absence of the white-noise limit. Some of 
the population growth models subjected to a white environmental noise changes the 
population-size dependence of the mean time to extinction (MTE) from an exponential 
to a power-law with a large exponent [5]. The introduction of the colored Gaussian noise 
changes this exponent, reducing it at a hxed noise magnitude. For a long correlation 
time of the environmental noise the MTE becomes independent of the population size 
for a strong enough noise [6]. 

Investigation of the effects caused by the presence of a colored noise in physical 
systems has some practical implications. Study of colored-noise-induced synchronization 
in chaotic systems indicates that the critical amplitude required for synchronization is 
generally smaller for the white noise as compared with the colored noise [7] . A practical 
implication is that, in situations where synchronization is undesirable, a simple control 
strategy is to place hlters in the system so as to make the noise source as colored 
as possible. In the systems exhibiting the phenomenon of stochastic resonance an 
exponentially correlated noise (“red” noise) leads to a reduction of signal amplihcation 
and the peak of stochastic resonance moves to a larger noise intensity when the 
correlation time increases |8]. “Pink noise” or 1// noise, as opposed to white noise also 
leads to a reduction of signal amplihcation, but resonance peak arises for lower noise 
intensity, if special conditions are satished [9] . This is important for understanding weak 
signal transmission trough noisy environments. 

Signals having the PSD at low frequencies / of the form S{f) 1//^ with /? 
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close to 1 are commonly referred to as “pink noise” or “1// noise”. Such signals are 
often found in physics and in many other helds uni [m ini tta El Ea Eg. Since the 
discovery of 1// noise numerous models and theories have been proposed, for a recent 
review see nzi. Mostly 1// noise is considered as a Gaussian process [HI EH], but 
sometimes the 1// fluctuations are non-Gaussian [201 El]- The Brownian motion as a 
source of 1// noise was hrst proposed in the seminal paper by Marinari et al. [22], 
where it was suggested that 1// noise can result from a random walk in a random 
environment. Starting from the model of 1// noise as a Brownian motion of inter¬ 
pulse durations [23l El [25] , nonlinear SDEs generating signals with 1// spectrum have 
been derived in [261127]. A special case of this nonlinear SDE has been obtained using 
Kirman’s agent model [28]. Such nonlinear SDEs have been used to describe signals in 
socio-economical systems [2^ 1^ . In this paper we consider the motion of a Brownian 
particle in an inhomogeneous environment and described by the nonlinear SDE yielding 
1/f spectrum. 

The influence of the colored noise on the dynamics of a Brownian particle 
immersed in a fluid where a temperature gradient is present can lead to interesting 
phenomena. The particle can exhibit a directed motion in response to the temperature 
gradient. Furthermore, study of stationary particle distribution shows that particles 
can accumulate towards the colder (positive thermophoresis) or the hotter (negative 
thermophoresis) regions depending on their physical parameters and, in particular, on 
the dependence of their mobility on the temperature [31]. The velocity of this motion 
can vary both in magnitude and sign, as observed in experiments [32]. However, in 
this case, no external force is actually acting on the particles [33] • Theoretical models 
suggest [31] the presence of a colored noise, as opposed to a white noise, is crucial for the 
emergence of such thermophoretic effects. Analysis of the steady-state dynamics of an 
overdamped classical particle in asymmetric multidimensional potential driven by the 
noise with an arbitrary correlation function has shown that the correlated noise breaks 
detailed balance, thereby exploiting the spatial asymmetries in potential to produce 
local drifts and rotations [3l]. These interesting Endings motivated us to investigate 
the motion of a Brownian particle in an inhomogeneous environment and subjected to 
a colored noise, as opposed to a white noise. 

The paper is organized as follows: In section [2] we review nonlinear SDEs driven by 
a white noise and yielding power-law steady state probability density function (PDF) of 
the generated signal. We estimate when the signal generated by such an SDE has 1/f 
PSD in a wide region of frequencies. In section[3]we show that nonlinear SDEs generating 
power-law distributed precesses with 1//^ spectrum can result from difiusive particle 
motion in inhomogeneous medium. In section E] we present methods that we use to 
study the influence of colored noise and in section [5] we investigate the effect of colored 
noise in our model: numerically solve obtained equations and compare the PDF and 
PSD of the signal with analytical estimations. Section [6] summarizes our Endings. 
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2. Nonlinear stochastic differential equations generating signals with 1// 
spectrum 


The nonlinear SDEs generating power-law distributed precesses with 1//^ noise have 
been derived in papers [26l [27l |35] starting from the point process model [23l EH 1 ^ . 
The general expression for the proposed class of ltd SDEs is 

dx = a^ L- j x^^-^dt + ax^dWt. ( 1 ) 


Here x is the signal, ?7 7 ^ 1 is the exponent of the power-law multiplicative noise, A 
defines the exponent of the power-law steady-state PDF of the signal, Wt is a standard 
Wiener process (the Brownian motion) and a is a scaling constant determining the 
intensity of the noise. The Fokker-Planck equation corresponding to SDE ([1]) gives the 
power-law probability density function (PDF) of the signal intensity Po{x) ~ x~^ with 
the exponent A. The non-linear SDE ([1]) has the simplest form of the multiplicative 
noise term, ax'^dWt- In papers [28l[36] the nonlinear SDE of type ([T]) has been obtained 
starting from a simple agent-based model describing the herding behavior. 

ltd SDEs are typically used in economics [37] and biology [38|. On the other 
hand, Stratonovich SDEs are more suitable for real systems with correlated, non¬ 
white noise, for example, for noise-driven electrical circuits [39] • The Stratonovich SDE 
corresponding to ltd SDE ([1]) is [40] 

dx = “ A)x^^“^dt -I- ax^ o dWt . (2) 

Note, that the choice of Stratonovich or ltd convention depends not only on the 
correlation time of the noise but also on the delay in the feedback [41] . 

For A > 1 the distribution Po{x) diverges as x —)• 0, therefore the diffusion of 
the stochastic variable x should be restricted at least from the side of small values, or 
equation ([1]) should be modihed. The simplest choice of the restriction is the reflective 
boundary conditions at x = Xmin and x = Xmax- Another choice would be modification 
of equation ([1]) to get rapidly decreasing steady state PDF when the stochastic variable 
X acquires values outside of the interval [xmin, a^max]- For example, the steady state PDF 


Po{x) -rcxp 


x^ 


/y>' ^ 
O/ri 


X 


X'' 


x'‘ 


( 3 ) 


with m > 0 has a power-law form when Xmin -C x <C x^ax and exponential cut-offs 
when X is outside of the interval [xmin, a^max]- Such exponentially restricted diffusion is 
generated by the SDE 


dx = 


'2 2 V I" 


X 


r^Ti 


x^a-Mt + ax^dWt 


( 4 ) 


obtained from equation ([T]) by introducing the additional terms in the drift. 

The PSD of the signals generated by the SDE ([1]) can be estimated using the 
(approximate) scaling properties of the signals, as it is done in [42]. Since the Wiener 
process has the scaling property dWat = changing the variable x in equation 
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Figure 1. (a) The steady-state PDF of the signal generated by equation © with 
reflective boundaries at Xmin and ccmax- The dashed line shows the power-law with the 
exponent —3. (b) The PSD of such a signal. The gray line shows the slope f~^. Used 
parameters are ry = 2, A = 3, Xmin = !> a^max = 1000, and (7 = 1. 


(II]) to the scaled variable Xs = ax or introducing the scaled time tg = one gets 

the same resulting equation. Thus change of the scale of the variable x and change 
of time scale are equivalent, leading to the following scaling property of the transition 
probability (the conditional probability that at time t the signal has value x' with the 
condition that at time t = 0 the signal had the value x): 

aP{ax', t\ax, 0) = P{x', a^t\x, 0), (5) 

with the exponent /i being ja = 2{r] — 1). As has been shown in [12], the power-law 
steady state PDF Po{x) ~ x~^ and the scaling property of the transition probability ([5|) 
lead to the power-law behavior of the PSD 

in a wide range of frequencies. 

The presence of the restrictions of diffusion at x = Xmin and x = Xmax makes the 
scaling ([5|) not exact and this limits the power-law part of the PSD to a hnite range of 
frequencies fmin <^ / <^ /max- The frequency range where the PSD has 1//^ behavior 
is estimated in |12] as 

< 27r/ < ,//>!, (7) 

< 27r/ < CX^X , 7^ < 1 . 

This equation shows that the frequency range grows with increasing of the exponent 
T], the frequency range becomes zero when rj = 1. By increasing the ratio Xmax/a^min 
one can get arbitrarily wide range of the frequencies where the PSD has 1//^ behavior. 
Note, that pure 1//^ PSD is physically impossible because the total power would be 
inhnite. Therefore, we consider signals with PSD having 1//^ behavior only in some 
wide intermediate region of frequencies, fmin <^ / /max, whereas for small frequencies 
/ /min the PSD is bounded. 
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Figure 2. The schematic representation of an inhomogeneous medium. The non¬ 
equilibrium system is described as a large number N of regions at local equilibrium, 
each region having different inverse temperature j3i. When the number of regions is 
sufficient large we can use approximation (1281) . 

For A = 3 we get that (3 = 1 and SDE ([T]) gives signal exhibiting 1// noise. 
Comparison of the nnmerically obtained steady state PDF and the PSD with analytical 
expressions for SDE ([T]) with rj = 2 and A = 3 is presented in hgnre [TJ For 
the nnmerical solntion we nse the Euler-Marnyama approximation, transforming the 
differential eqnations to difference eqnations. We nse a variable time step, decreasing 
with the increase of x. As in |26l ET] we choose the time step in snch a way that 
the coefficient before noise becomes proportional to the hrst power of x. Very similar 
nnmerical resnlts one gets also by using the Milstein approximation |35] . We see a good 
agreement of the numerical results with analytical expressions. A numerical solution of 
the equations conhrms the presence of the frequency region for which the power spectral 
density has 1//^ dependence. The 1// interval in the PSD in hgnre [1] is approximately 
between /min ~ 10° and /max ~ 10° and is much narrower than the width of the region 
1 -C / -C 10° predicted by equation ([7]). The width of this region can be increased by 
increasing the ratio between the minimum and the maximum values of the stochastic 
variable x. In addition, the region in the power spectral density with the power-law 
behavior depends on the exponent r/: if r/ = 1 then this width is zero; the width 
increases with increasing the difference — 1| |43j . 

3. Nonlinear SDE resulting from motion in inhomogeneous medium 

The nonlinear SDE ([T]) generating signals with 1// spectrum in a wide range of 
frequencies has been used so far to describe socio-economical systems |29l |30]. The 
derivation of the equations has been quite abstract and physical interpretation of 
assumptions made in the derivation is not very clear. In this Section we present a 
physical model where such equation can be relevant. We expect that this derivation 
leads to a better understanding which systems can be described using equation ([T]) or 
equation (|2]). 
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We will consider the Brownian motion of a small macroscopic particle in 
an inhomogeneons medinm. We assnme that this medinm has reached local 
thermodynamical eqnilibrinm bnt not the global one and the temperatnre can be 
considered as a fnnction of coordinate. Schematically snch a medinm is shown in 
hgnre [2l Brownian motion of small macroscopic particles in a liqnid or a gas results 
from unbalanced bombardments due to surrounding molecules. Usually the Brownian 
motion is described by a Langevin equation that includes the influence of the “bath” of 
surrounding molecules as a time-dependent stochastic force that is commonly assumed 
to be a white Gaussian noise. This assumption is is valid when the correlation time 
of fluctuations is short, much shorter than the time scale of the macroscopic motion, 
and the interactions with the bath are weak. The effect of large correlation time of 
fluctuations will be considered in next Section. In the case of strong collisions between 
the particle and the surrounding environment the noise is non-Gaussian and we have 
so called Levy flights [H]. Nonlinear SDEs Levy noise and generating signals with 1// 
spectrum are proposed in ^5] . 

Langevin equations for one-dimensional motion of a Brownian particle are |16] 




(9) 


Here x is the coordinate and v is the velocity of the Brownian particle, m is the mass 
of the particle, 7 is the relaxation rate and ^{t) is the (5-correlated white noise. In 
general, equations similar to (jS]), ([ 2 ]) can be used to describe a variety of systems: noisy 
electronic circuits, laser light intensity fluctuations |16] and others. We assume that 
there is temperature gradient in the medium, therefore the inverse temperature /5(a;) 
depends on coordinate x. In the case when (3{x) = (3 = const, equations (jH]) and ([9]) 
describe Brownian motion in the medium with constant temperature T = j where 
k-Q is Boltzmann constant. 

In high friction (also called overdamped) limit [iQ] the relaxation rate is large, 
7 S> \dv/dt\. Performing adiabatic elimination of the velocity as in |17], we obtain the 
equation 



( 10 ) 


Here (3'{x) = d(3{x)/(ix. This SDE should be interpreted according to the Stratonovich 
convention. Note, that the second term in the right hand side of equation (ITOH arises 


due to position dependence of the stochastic force in equation ([H]) [47]. The Ito SDE 
corresponding to flTOl) is 



( 11 ) 
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For calculating stationary distribution of position x in high friction limit we will use a 
Fokker-Planck equation instead of SDE ffTOj) . The Fokker-Planck equation corresponding 
to (na can be written as 




d 19^ P{x, t) 

-F{x)P{x,t)+ ^ ^ 


'ymdx ^ ' 'jmdx‘^ /3{x) ^ ^ 

By setting the time derivative to zero we obtain an analytical expression for steady state 
distribution Po{x)\ 


Po{x) = C/3{x)exp yj F{x')/3{x')dx'j . (13) 

Let us consider the situation when the dependence of the inverse temperature on 
the coordinate is described by a power law, 

^{x) = bx~^ . (14) 

Here 0 is a power law exponent and 6 is a constant. This is quite reasonable assumption; 
for example, if 6^ = 1 equation flTT)) represents a case where we have steady state 
heat transfer due to temperature difference T 2 — Ti between the beginning and the 
end of the system. We assume that there are reflective boundaries at x^m and a^max 
and the motion is limited to values of x between Xmin and Xmax- When 9 = 1 then 
the temperatures should obey the relation T 2 /T 1 = Xmax/a^min and the coefficient b is 
b = (xmax — Xrain)/kBiT 2 — Ti). This case is presented in £gure|2l 

The external force affecting the particle F{x) we express via the gradient of the 
potential V{x): 


F(x) = - — V(x). 
dx 

We choose the subharmonic potential proportional to the temperature: 


(15) 


(16) 


e ~J i3(x) 

For convenience we write the coefficient of proportionality as vjO — 1. As we will see 
in equation ffl^ . the parameter v gives the power law exponent in the steady state 
distribution of x. Taking into account equation (ITT)) we see that the potential has the 
power law form with the same exponent 9. The expression for the external force then is 

F(x) = ^x®-^ 

Using inverse temperature ca and force 03 the equation (na for the particle 
coordinate becomes 

d ... 9 (1 

2 


-x(f) =-- 

at 'ymb 


v 

— w I X 




+ X2 


-«(i) 


(18) 


9 J V T’^b 

By using equations ffT3|) . flTT)) and ffT7)l we obtain distribution of particles in high friction 
limit 


Pn(x) = 


z/- 1 


‘^min ‘^max 


X 


(19) 
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Figure 3. The steady-state PDF of the signal generated by the Langevin equations 
(HID, (H3D with reflective boundaries at Xj^in and imax and the parameters: (a) 0 = 1, 
V = 2), iinin = 1, imax = 1000, {b) 0 = 0, V = 3, imin = 1, imax = 100. The dashed line 
shows the power-law with the exponent —3. (c) The PSD of the signal corresponding 
to the parameters in (b) case. The gray line shows the slope f~^. 


Calculating distribution of particles we assumed that there are reflective boundaries at 
X min and Xmax and the motion is limited to values of x between x^m and Xmax- Equation 
(lT 8 l) has the same form as Stratonovich SDE ([21). 

3.1. Numerical solution 

To check the validity of the adiabatic elimination of the velocity, we solve the Langevin 
equations ([H]), (Ej) numerically. For numerical solution it is convenient to introduce 
dimensionless time t, position x and velocity v. When the inverse temperature /3{x) and 
the force F{x) are given by equations (1TT|) and (fT7|) . we can use i = 'jt and 


X = 2-ex , 

(20) 

V = 7 “^( 7 ^m 6 ) 2 ^n . 

(21) 

The equations ((HD, (ED in the dimensionless variables have the form 


-in = —V + {9 — iy)x^~^ + V2x^^{i ), 

dt 

(22) 

d 

—X = V . 

dt 

(23) 

The reflective boundaries at Xmin and Xmax become 


^min (T 

(24) 

^max (T • 

(25) 

The requirement of large friction 7 , necessary for the adiabatic 
velocity, leads to the requirement Xminj^max 3> 1 when 9 <2. 

elimination of the 

Applying the Euler scheme with the step h to (j22il. (l23D yields the following 

equations: 


Vk+i = Vk - Vkh + {9- n)xl~^h + sj2xlh^k , 

(26) 

Xk+l = Xk + Vkh . 

(27) 


Here are independent Gaussian random variables with zero mean and unit variance. 
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As an example we solve the Langevin equations 0221) . fl23p with reflective boundaries 
at Xjnirn ^max and the parameters n = 3, 9 = 1 or 9 = 0. The value of the parameter 
9 = 0 means that the temperature is constant. The steady state PDF of the particle 
position Po{x) and the power spectral density S{f) are presented in figure [3l In figure [3^ 
we see a good agreement with the distribution of particles in high friction limit flT^ . 
This confirms the validity of the adiabatic elimination of the velocity. Figure [3l3 confirms 
the presence of the frequency region with 1// behavior of the power spectral density, 
the 1// interval in the PSD is approximately between /mm ~ 10“'^ and /max ~ 10“^. 
The width of this interval can be increased by increasing the ratio Xmax/^min- The width 
of 1//^ region in the PSD also increases with increasing of 1^/2 — 1|. 

Not only coordinate fluctuations yield power law PSD. The Langevin equation with 
postion dependent temperature (|T4|) and external force (|T7jl can also lead to the power 
law PSD of the absolute velocity fluctuations. 


3.2. Superstatistics and velocity fluctuations 


In the motion of the particle there are two time scales: the scale in which the particle 
is able to reach a local equilibrium and the scale at which the fluctuating temperature 
changes. We assume that temperature fluctuations are slow, that is, the time scale 
at which observable change of temperature happens is much larger that relaxation 
time of the particle I/7. This assumption empowers us to derive simpler equations 
for particle velocity and apply the superstaticlical approach |1HI SHI EHl EU [52] to obtain 
statistical properties of particle velocity fluctuations. The superstatistical framework 
has been successfully applied to a wide range of problems, such as interactions between 
hadrons from cosmic rays [53|, fluid turbulence [SU [55l EH EZj, granular material [58] 
and electronics [59] • In the long-term the nonequilibrium system is described by the 
superposition of different local dynamics at different time intervals. Superstatistics 
is a description of the complex system under consideration by a superposition of two 
statistics, one corresponding to ordinary statistical mechanics (Langevin equation) and 
the other one corresponding to a slowly varying parameter, in this case the inverse 
temperature fl{x). In the superstatistical approach the distribution of the velocity of 
the particle is 

^ /•/3max 

P{v) = J2fWPWi)^ fW)P{vW)d(l, (28) 

i=\ ^/^min 

where P{v\fl) is a local velocity distribution for the particle near the position x when the 
temperature can be assumed to be constant and equal to fl{x) = (3. For the Langevin 
equation in the form of equation ([H]) the local velocity distribution is 


P{v\fl) 



(29) 


The function f{fl) is the distribution of the inverse temperature. The distribution of 
the inverse temperature f{fl) can be found from the stationary distribution Pq{x) of 
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coordinate x by using the relation [50] 

f(P) = (^) ' PoW))- 

The parameters 

Mmax ’ ymm 


(30) 

(31) 


are maximal and minimal inverse temperatures. From equations 
obtain the distribution of inverse temperature 

n-1 


m = 


e 


Rn 


■f3- 




■ (IT9i) and (ITD we 
(32) 


Now we will describe fluctuations of particle velocity. For this purpose we introduce 
the local average of absolute value of the velocity 


v= \v\P{v\(3)dv. 


(33) 


Keeping the assumption that local distribution of the velocity remains Gaussian fl29|) . 
we have 

2 

7Tml3{x) 

By using equation flMj) and changing the variable from a; to h in Stratonovich SDE flTSl) 
|4Uj . we obtain an SDE for the average velocity 

= ^a^{ri - . (35) 

at 2 

Here we introduced the new parameters 

^ = + l (36) 

and noise intensity 
a = 6 

The equation fl35l) for the average absolute velocity is identical to SDE ([2]) generating 
power-law distributed signals with l//“ spectrum. The multiplicity of noise rj in 
equation (|35ll depends only on power-law exponent 9, whereas the exponent A of the 
power law part of distribution of v depends on both parameters 6 and u. According to 
equation ([6]) we have 1// noise when X = 3oYn = 9 + l. One can expect that fast 
velocity fluctuations do not influence the spectrum at small frequencies and 1// noise 
is visible not only in the spectrum of the local average v but also in the spectrum of the 
absolute value of the velocity |n|. Numerical calculation conhrms this expectation. 



(37) 
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Figure 4. (a) The steady-state PDF of the absolute velocity obtained by numerically 
solving the Langevin equations (l2^ . (l23l) with reflective boundaries at imin and imax- 
The gray line shows the PDF from equation (1381) multiplied by 2. (b) The PSD of the 
absolute velocity fluctuations. The gray line shows the slope The parameters 

used are: 6 = 1, i/ = 2, imin = Ij imax = 10^. 


C(t>) = 


From equations 
A-1 


and 


we get 


2^/F 


A-1 

max 


A 


A-1 

2 

min 


A-1 

2 


l-A 


A 1 2. 

2 ’ 2 ^^ 


r ( 2 ’ 2 Amax 


(38) 


Here r(a, z) = is the incomplete gamma function. When -C 1 

and ^mn^/dmax 3> 1 then from equation (j38|l it follows that the distribution of velocities 
has a power law form P{v) ~ 1^1”^- 

Thus we obtain 1// noise in the fluctuations of the absolute value of the velocity 
when the velocity distribution has a power-law part P{v) ~ v~^ and temperature 
distribution is flat, f{/3) = const. When 6 = 1 we have a simple system exhibiting 1// 
fluctuations: the system consists of a Brownian particle affected by a linear potential 
V (x) and moving in the medium where steady state heat transfer is present due to the 
difference of temperatures at the ends of the medium. 

As an example we solve the Langevin equations fl2^ . fl23|) with reflective boundaries 
at £min, almax and the parameters n = 2, 9 = 1. The steady state PDF of the absolute 
velocity Ddhl) and the power spectral density S{f) of the absolute velocity fluctuations 
are presented in hgure IH The steady state PDF Ddnl) is twice as large as the PDF 
given by equation because both positive and negative velocities with the same 
absolute value occur with equal probabilities. The PDF of the dimensionless velocity v 
can be obtained by setting m = 1, 6 = 1. In hgure 0^ we see a good agreement with 
the analytical expression for the steady state PDF. Figure |Dd conhrms the presence of 
the frequency region with the power law behavior of the PSD. The power law region 
in the PSD is approximately between ~ 10“^ and /max ~ 10“^. However, there is 
a slight difference between the power law exponent equal to 1 predicted by analytical 
expressions and the numerically obtained power law exponent 1.1. This disagreement 
arises because the analytical prediction of the power law exponent in the spectrum is 
only approximate. 
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= f{r)+3(x)e(t). 


(39) 


A system subjected to the colored noise is described by the Langevin equation with a 
time-dependent stochastic force that includes the influence of the “bath” of surrounding 
molecules: 

dx 

dt 

It is a well known result that if we approximate white noise by a smooth, colored 
process, then at the limit as the correlation time of the approximation tends to zero, the 
smoothed stochastic integral converges to the Stratonovich stochastic integral [61]. We 
assume that the stochastic force in equation fl3^ is a Gaussian noise having correlation 
time comparable with time scale of the macroscopic motion and that e{t) is Markovian 
process. In such a case the noise e satishes Ornstein-Uhlenbeck process with exponential 
correlation function [62|. Thus we describe the system perturbed by a colored noise as 
two-dimensional Markovian flow: 


^ = f{x)+g{x)e, 


de 

dt 


= —e ■ 

T 




T 


m- 


(40) 

(41) 


Here ^{t) is the white noise, {^it)^{s)) = 6{t — s), the parameter r is the correlation 
time of the colored noise and D is the noise intensity. The autocorrelation function of 
the colored noise is 


{e{t)e{s)) = — exp 

T 




T 


(42) 


It is possible to write two dimensional Fokker-Planck equation for equations (140]) . 
flTT]) and obtain two dimensional P{x, e) density as its solution. However, usually it is 
enough to know the distribution P{x) of the signal x, which can be formally obtained 
by averaging P{x,e) over the noise e. It is problematic even get approximate analytical 
solutions for P{x,e) [63]. More convenient way is to get P{x) from an approximate 
Fokker-Planck equation just for one variable. Such equation can be obtained by using 
the unihed colored noise approximation [M] , which is a form of adiabatic approximation. 
To make this approximation we eliminate the variable £ and get the equation 

,2. ^'(X) 


d X 


dx\ /I .g'{x)\ dx I j., . 


+ —g(x)((t) ■ (43) 


dP g{x) 

We assume that the variable x changes slowly and drop small terms containing d^x/dP 
and (dx/dt)^, obtaining the approximate equation 

dx /(x) y/2Dg{x) 




1-r (/'(x) -/(x)^) 


m . (44) 


l-r(^/'(x)-/(x)^ 

This equation should be interpreted in the Stratonovich sense. In more general case, 
when we cannot neglect inertia and drop the second derivative di^x/dP, the question 
whether the equation obtained using adiabatic approximation should be interpreted 
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in Ito or Stratonovich sense still remains an open question, so called Ito-Stratonovich 
problem |65]. However, at least for specific systems in white noise limit, it can be 
determined which interpretation is correct. For example, it has been shown for a 
simplihed model of the preferential concentration of inertial particles in a turbulent 
velocity held [66], that the equation obtained using adiabatic elimination in white noise 
limit became the Stratonovich equation [67]. The Stratonovich interpretation should 
be used if the correlation time of the noise is much larger than the relaxation rate of 
the system. In an opposite case the equation should be interpreted in ltd sense. If 
relaxation rates are of the similar magnitude as the correlation time, we get an equation 
with noise induced drift that is different from Stratonovich drift. 

The Fokker-Planck equation corresponding to the Stratonovich equation dx/dt = 
fc{x) +gc{x)^{t) is [40] 

^P(a:, t) = —^fc{x)P{x, t) + ^-^g,{x)-^gc{x)P{x, t ). (45) 

The applicability of equations fl44|) and fl4^ has limitation due to neglect of higher order 
derivatives in equation (1431) . These equations describe the dynamics correctly [6l] for 
times t obeying 

T 


t > 


1-T { fix) - f{x)^ 


and in the space regions obeying 

V2DTg{x) 


f'{x) 




9{^) 


fix) 


fix) 


< 1 . 


(46) 


(47) 


5. Influence of colored noise on the stochastic differential equation 
generating signals with 1// spectrum 


If the nonlinear SDE generating signals with 1// spectrum is a result of a Brownian 
motion in an inhomogeneous medium then the finite correlation time of the “bath” can 
become important. In this section we use the results presented in section |3]to investigate 
the influence of the colored noise. Instead of white noise we add colored noise e(t) to 
the Stratonovich equation ([2]) obtaining the equations 


^ = ^fiv - A)x^'' ^ + ax'^e{t ), 

de 1 X 

= —e + -^(t). 


dt T T 
After unified colored noise approximation 

dx _ |a^(r7 - A)x^^“^ 

df 1 


we get 


+ 


ax ‘ 


\Ta 


— l){rj — X)x‘^'x-d 1 


\Ta 


‘^{jj — l){ri — X)x‘^^'^~d 


W)- 


(48) 

(49) 

(50) 


If T is large then equation (l50|) has a simpler form 

dx X 2 x^~^ 

dt rin — 1) rain — 1)(A — rt) 


m- 


fg-l) Ta{r]-l){X-g) 


(51) 
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Equation fl5D]) should be interpreted in the Stratonovich sense. Converting to ltd 
interpretation [iQ] we have 


dx = 


1 

2 7(x) 


^ , 2-7 ^ 2(p-l) 
V ~ —TT 

7(x) 7(x)^ 


dt + 


ax' 

7(x) 


dlE* 


where 


7(x) = 1 — -Ta^{r] — l){r] — 


According to equation fH7|) . approximation fl50|l is valid when 
^/Ta\2r] — 


1 — ^Ta'^iji — l){ri — 


< 1 


Steady state PDF corresponding to equation (jHIil) is 


Pq{x) ~ X 


-A 


1 2. 
-ra (ri 
2 ^' 


1){t] — \)x‘^^^ 


exp 


1 2. 
iv 


X)^x 


2^2(r,-l) 


(52) 


(53) 


(54) 


(55) 


1 2 | 
-ra Iri 
2 ' ' 


We see that the colored noise introduces an exponential cut-off in the steady state 
PDF Po{x) and naturally limits the range of diffusion of the stochastic variable x. The 
exponential cut-off is at large values of x when 7 > 1 and at small values of x when 
7 < 1. 

Comparing equation (1501) with ([2]) we see that the influence of the hnite correlation 
time r of the noise can be neglected when 

1||7 — -C 1. (56) 

Let us consider the case 7 > 1. Then according to equation fl56|) . the influence of the 
hnite correlation time r of the noise can be neglected when x <C Xr, where 

2 1 2 (> 7 - 1 ) 

^ [ra2(7-l)|7-A|J 

If the diffusion is restricted to the region Xmin < x < x^ax then the spectrum has a 
power-law part in the frequency range given by ([7]). If x,- > Xmax we expect no change 
in the power-law part of the spectrum. If Xr < Xmax then, replacing the maximum value 
of X by Xt we get that the replacement of the white noise by the colored noise leaves 
the power-law part of the spectrum in the frequency range 


2 2(»;-l) 
^ ^min 


< 27r/ -C 


(58) 


r(r(-l)|i;-A| 

If 7 < 1 then the inhuence of the hnite correlation time r of the noise can be 
neglected when x 3> x,-, where Xr is given by flFT)) . When the diffusion is restricted to 
the region x^in < x < Xmax, we expect no change in the power-law part of the spectrum 
when Xr < Xmin- If Xr > Xmin then replacing the minimum value of x by Xr in equation 
(]7|) we can estimate that the power-law part of the spectrum should be in the frequency 
range 

2 


2 - 2 ( 1 -,?) 
^ ‘^max 


< 2'Kf -C 


(59) 


t( 1 - )))|r)- A| ' 

Thus the introduction of the colored noise into equation (j2]) can narrow the range of 
frequencies where the PSD behaves as 1//^ by decreasing the upper limiting frequency. 
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t t 

Figure 5. (a) Typical signal generated by equations (H51) . (H^l) with r = 0.01. (b) 
Typical signal generated by equation @ corresponding to r = 0 (white noise). We used 
reflective boundaries at Xmin = 1 and Xmax = 1000. Other parameters of equations are 
77 = 2, A = 3, and ct = 1. 






Figure 6. (a,c) The steady-state PDF of the signal generated by equations (IS51) 1- 
(EZl) with reflective boundaries at Xmin and Xmax- The dashed line shows analytical 
approximation (1551) . (b,d) The PSD of such a signal. The gray line shows the PSD 
of the signal generated by equation The correlation time is (a,b) r = 10 (c,d) 

r = 10“^. Other parameters are 77 = 2, A = 3, Xmin = 1, a^max = 1000, and cr = 1. 


5.1. Numerical solution 

To check the validity of the approximations we performed numerical solution of equations 
(HE]), For numerical solution we used the Euler scheme |6Hll69]. Applying the Euler 
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scheme with the step h to 


yields the following equations: 


^k+l “1“ 


1 2. 


X)xl'^ ^ + axlzk ) h , 


Zk+l — ^k 


Izkh + 

T T 


(60) 

(61) 


Here are independent Gaussian random variables with zero mean and unit variance. 
For rj < 1 and for negative rj it is convenient to use the Euler scheme with a constant 
time step h. However, in the case of ?7 > 1 at large values of x the coefficient in equation 
(IHnjl becomes large and thus requires a very small time step. A more effective solution 
is to use a variable time step, decreasing with the increase of x, as has been done in 
Refs. [26l[27]. Variable time step hk = nT/{ax^) leads to the equations 

'1 


Xk+i =Xk + KT { -a{r] 


X)xl ^ + Zk 


^k+l ^k 


K 


tk+1 — tk X- 


KT 


fZk + 


K 


arx 


^k y 


(TXl 


(62) 

(63) 

(64) 


Here k -C 1 is a small parameter. 

As an example, we solve equations fH8l) . with the parameters rj = 2, X = 3, 


a = 1 and reflective boundaries at x^ 
get 


= 1 , Xr 


= 1000. From equations 


Xk+l = Xk + KT 


^k 2 


K 


Zk+1 — ^k 


tk+l — + 


KT 


rZk + 


Xk 


^k y 


Xf 


we 

(65) 

( 66 ) 

(67) 


The generated signal is shown in hgure [5l We see that the hnite correlation time r of 
the noise leads to a smoother signal compared to the equation with r = 0. The steady 
state PDF Po{x) and the power spectral density S{f) for two different values of r are 
presented in hgure [6l From hgure |6^ we see that the unihed colored noise approximation 
correctly predicts the exponential cut-oh in the steady state PDF at large values of x, 
although the actual position of the cut-oh slightly dihers from the cut-oh predicted by 
equation (|55il . As hgure |6b shows, the presence of the hnite correlation time r makes 
the power-law part in the spectrum narrower. The upper limiting frequency of the 
power-law region grows with decreasing of r, as is qualitatively predicted by equation 
(1581) . The steady state PDF and the PSD of the generated signal corresponding to much 
smaller value of the correlation time r are shown in hgure [6)D,d. For this value of r the 
exponential cut-oh due to hnite correlation time is larger than the upper boundary a:max, 
thus we see almost no diherences from the case of uncorrelated noise, r = 0. 
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6. Discussion 

The nonlinear SDE ([T]) generating signals with 1// spectrum in a wide range of 
frequencies has been used so far to describe socio-economical systems |29l |30]. The 
derivation of the equations has been quite abstract and physical interpretation of 
assumptions made in the derivation is not very clear. In this paper we propose a physical 
model where such equations can be relevant. This model is described by Stratonovich 
(E]) instead of ltd ([I]) SDE and provides insights which physical systems can be described 
by such nonlinear SDEs. 

We have shown that nonlinear SDEs generating power-law distributed precesses 
with 1//^ spectrum can result from diffusive particle motion in inhomogeneous medium. 
The SDE (l35|) for velocity fluctuations or SDE (fT8|) for particle coordinate are simplihed 
versions of Langevin equations (IH]), (E]) for one-dimensional motion of a Brownian 
particle. We neglected viscosity dependence on temperature and inertia of particle. 
In general, equations similar to these can be used to describe a variety of systems: 
noisy electronic circuits, laser light intensity fluctuations |16] and others. We assumed 
that the inverse temperature (3{x) depends on coordinate x and this dependence is 
of a power law form. Such a description is valid for a medium that has reached 
local thermodynamical equilibrium but not the global one, and the temperature can 
be considered as a function of coordinate. The power law dependence of the inverse 
temperature (3{x) on the stochastic variable x can be caused by non-homogeneity of a 
bath. This non-homogeneity can arise from a complex scale free structure of the bath 
as is in the case of porous media uni or from the bath not being in an equilibrium. 

In high friction limit, if the particle is affected by a subharmonic potential 
proportional to the local temperature, the motion of the particle can be described 
by the equation similar to equation ([2]). For example, we can consider a Brownian 
particle affected by a linear potential V{x) and moving in the medium where steady 
state heat transfer is present due to the difference of temperatures at the ends of the 
medium. From the properties of equation (j2]), presented in section Ej it follows that 
the spectrum of the fluctuations of the particle position x(t) in such a system can 
have a frequency region where the spectrum has a power-law behavior. The width 
of this frequency region increases with the increase of the length of the medium in 
which the particle moves. Equation ([2]) can also describe the fluctuations of the local 
average of the absolute value of the velocity, if temperature fluctuations are slow and 
the superstatistical approach can be used. We obtain 1// noise in the fluctuations of 
the absolute value of the velocity when the velocity distribution has a power-law part 
P{v) V ^ and temperature distribution is flat, f{(3) = const. 

The correlation of collisions between the Brownian particle and the surrounding 
molecules can lead to the situation where the hnite correlation time becomes important, 
thus we have investigated the effect of colored noise in our model. Using the unihed 
colored noise approximation we get that the hnite correlation time leads to the additional 
restriction of the diffusion. Existence of colored noise leads to an exponential cut-off of 
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the PDF of particle positions either from large values when rj > 1, or from small values 
when rj < 1. Such a restriction of the diffusion is a result of the multiplicative colored 
noise in equation fHSll . Narrower power law part in the PDF of the particle positions 
results in the narrower range of frequencies where the spectrum has 1//^ behavior. 
When rj > 1 , the end of the power-law part of the spectrum at large frequencies is 
inversely proportional to the correlation time r of the noise. However, for sufficiently 
small correlation time, when the restriction of the diffusion due to colored noise is 
larger than the upper boundary a^max of the medium, the effects of the colored noise are 
negligible (see hgure [6]) and the properties of the signal do not differ signihcantly from 
the white noise case. 
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